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Abstract 

Two flavours Skyrmion will be analyzed analytically, in case of static and rotational 
Skyrme equations. Numerical solution of a nonlinear scalar field equation, i.e. the 
Skyrme equation, will be worked with finite difference method. This article is a more 
comprehensive version of SU(2) Skyrme Model for Hadron which have been published 
at Journal of Theoretical and Computational Studies, Volume 3 (2004) 0407. 
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1 Introduction 

Nonlinear physics is a general phenomenon in physics. It includes very long range from 
e.g. particle, nuclear, condensed matter, fluids, plasmas, biophysics (for example: nonlinear 
diffusion-reaction, DNA), to cosmology. 

Soliton is defined as classical solution of nonlinear wave equation which has properties: 
finite total energy, localized, nondispersive, stable with its profile of energy density distribution 
is like pulse centered in finite space pQ. 



This topics was presented at WCS2K5 UGM, Yogyakarta, Indonesia, 8 September 2005. 
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Skyrmion is a topological soliton in three dimensions. Following, we will discuss soliton as 
baryon model, especially Skyrme model for hadron which is formulated in extended nonlinear 
Sigma model [Tj, [2J. 



2 Topological Solitons 



Suppose that = is a differential equation or a system of differential equation involving a 
field or a set of fields denoted by <j) which are functions of D space coordinates x and a time 
coordinate t. We shall assume that = is the equation of motion which results from some 
Lagrangian filed theory. We also assume that the system has a corresponding "energy" 
and "energy density" E<f,{x,t) [3] 



where for all allowed field configurations <p, E^ is greater than or equal to zero. If is zero 
for all x, we call (ft a ground state or a vacuum solution, and denote it by (p vac . The vacuum 
solution need not be unique. 

Now suppose that <fi = <p c i is non-vacuum solution to F^ = 0. Following Coleman, we shall 
call 4> c i(x,t) a soliton solution if the following properties are satisfied [3]: 

(1) £0=^ is finite. 

(2) E^ = ^ cl {x,t) is nonsingular (finite) for all value of x and t, and localized for all times t. 
A solution is localized at any time t if there is a bounded region of space defined by 
E$=$ cl (x, t) > 6, 5 being any arbitrary number, fulfilling 



We say that the solution is localized for all times t if the bounded region can be chosen 
independent of t. 

(3) cf) c i(x,t) is nonsingular. 

(4) <f) c i(x,t) is nondissipative. 

Concerning (4), a solution 4> c i(x,t) is considered dissipative if 



In our discussions, we will be concerned primarily with "static" solutions for which we 
denote 4> c i{x,t) by $ d (£,£). 

// we exclude the vacuum, then static, nonsingular, localized solutions are automatically 
nondissipative. Thus for static solitons, we only require (1), (2) and (3). 

We will, however, be interested in static solutions which satisfy yet another requirements; 
namely that they be 




(1) 



< 5 < maxj E^ =4)c[ (x,t). 



(2) 



limt^ooHiax^ E^^Xjt) = 0. 



(3) 
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(5) classically stable. We call a static solution "classically stable" if under all possible 
infinitesimal fluctuations 5(f) of the fields about the static solution, 

£4> cl +s<t>>£<t> cV (4) 

Variational modes 5o0 which leave the energy unchanged are said to be "zero" frequency 
modes" or "zero modes". Thus, for zero modes, e^+s^ = e^ cl . Zero modes are generally 
associated with various symmetries of the soliton. Excluding the zero modes, a classically 
stable solution is a local minimum of the energy. 

Although the question of stability is ultimately a dynamical one, topology often plays an 
important role. We shall now elaborate on this point. Let Q be the set of all finite energy 
and nonsingular field configurations at some fixed time t. It is thus the configuration space 
of the system. A subset Qi of Q is said to be (path-) connected if any field <\>\ of Qi can be 
continuously deformed to any other field <fl x in Q x . 

Two subsets Qi and Q2 of Q are said to be disconnected if any field 0i of Qi cannot 
be continuously deformed to any field 02 in Q2- We shall consider the case where Q has 
N > 1 disconnected components Q n , each Q n being connected. Q is the union of all these 
disconnected components, 

Q = u^UQn- (5) 

Let 0(5?) and 4>'(x) be two fields in Q. <f)(x) and 4>'(x) are homotopic to each other, and 
write ~ 0', if there exists a sequence of fields <ffr\x), < r < 1, which is continuous 
in r and x, with <f>(°>(x) = 0(5f) and <fp-\x) = (f>'(x). With this definition, all fields within 
one component Q n are "homotopic" to each other, while a field belonging to Q n is not 
homotopic to a field 0(™'- ) belonging to Q n i when n' 7^ n [3]. 

We can treat all our examples assuming that n is countable. We make this assump- 
tion: What is the physical significance of this classification ? Consider the initial conditions 
(4>^ n \ d(j)( n > /dt) at time t = for the equations of motion. (These equations for simplicity are 
assumed to be second order in time). Assume that <p n ) G Q n . After a lapse of time T, suppose 
that (n) becomes /(n) and d<j) {n) /dt becomes d<p'^ n ' / dt. Since time evolution is assumed to be 
a continuous operation, it follows that <p'^ is homotopic to (j)( n \ Hence, (p 1 ^ G Q n . In other 
words, the value of n associated with the field cf)^ is a constant of the motion. Hence, if we 
define Qo to be the sector which contains the vacuum solution mc , then the configurations 
(j)( n > G Q n , 11 / 0, cannot be time evolved to mc or in fact, to any other 0( m ) G Q m with 
m ^ n. For such reasons, the sectors Q n ^o of field configurations are said to be topological 
stable. Soliton configurations which fulfill the five properties listed previously and which are 
in topologically stable sectors are called topological solitons [3]. 

3 Two Flavours Skyrmion as Hadron 

In particle physics, flavour is a quantum number of elementary particles. In quantum chro- 
modynamics (QCD), flavour is a global symmetry. In the electroweak theory, on the other 
hand, this symmetry is broken, and flavour changing processes exist [3]. 

If there are two or more particles which have identical interactions, then they may be 
interchanged without affecting the physics. Any (complex) linear combination of these two 
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particles give the same physics, as long as they are orthogonal or perpendicular to each other. 
In other words, the theory possesses symmetry transformations such as M Q), where u and d 
are the two fields, and M is any 2x2 unitary matrix with a unit determinant. Such matrices 
form a Lie group called SU(2), special unitary group. This is an example of flavour symmetry. 
The term "flavour" was first coined for use in the quark model of hadrons in 1968 [I]. 

In order to understand the charge radius of nucleon which have size roughly 1 Fermi, Tony 
Hilton Royle Skyrme in 1962 proposed idea that strongly interacting particles (hadrons) were 
locally concentrated static solution of extended nonlinear Sigma (Chiral) model. In (3+1)- 
dimensions of space-time, we will observe Skyrme model which describes hadrons as solitons 
(Skyrmions) from nonlinear Sigma (chiral) model field theory in internal symmetrical group 
of SU(2): two flavours Skyrmion. 

Skyrme's idea is unifying bosons and fermions in a common framework which provide 
a fundamental fields model consisted of only pion. The nucleon was obtained, as a certain 
classical configuration of the pion fields. Skyrmion is a topological soliton object, i.e. solution 
to the classical field equation with localized energy density. Various atoms should correspond 
to vortices of different connectivities in some underlying liquid. 

Following, some ideas which are proposed by Skyrme in relations with baryon model PQ, 

(1) Meson field can take its value in S3 manifold. As the result of this assumption, Skyrme 
discovers conserved quatity i.e. topological charge or winding number and he interpre- 
tates it as baryon number. 

(2) Solution of Skyrme field equation in spherical coordinate takes form: 

F a {r,0,4>) = g(r)n a {p,4>) (6) 

where g(r) is profile function which has spherical symmetric character and n a (a = 
1,2,3) is component of unit vector, n. F a (r,9,(p) is called Skyrme ansatz or hedgehog. 

(3) Skyrme ansatz describes stable extended particle with unit topological charge. 

(4) Skyrme ansatz with unit topological charge can be quantized as fermion, so it is possible 
to identify state of isotopic spin, /, and total spin, J, with nucleon doublet in case of 
I = J = ~, and with A resonance in case of / = J = |. 



4 SU(2) Skyrme Model 

SU(2) Skyrme model is very simple model, because it only consists of two flavours. This 
model is described by U — £/(x M ) function which has SU(2) group valued of (3+l)D space- 
time coordinates [TJ Q2J E] . 

The dynamics is determined by action 

S= [ d 4 x C, (7) 
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where 



C = Tr 



16 p ' 32a 2 
is the related Lagrangian density and 



^ V + J-j [L„ L„] [L M , LI + ^M^C/" 1 + 17-2/) 



l„ = u- l d u u 



(8) 



(9) 



are ie/t chiral currents, where F = 123 MeV is pion decay constant and a is dimensionless 
constant. The first term in equation |2j) is the SU(2) chiral model Lagrangian density, the 
second term is Skyrme term for stabilizing solitonic solution. The last term descibes the mass 
term where M n is pion (meson) mass. 

Euler-Lagrange equation of SU(2) Skyrme model is derived from least action principle 



5S = 0. 

If we take variation of action to equation (J^J), we obtain 

6S= f d 4 x SC. 



Substitute equation into (TTTT) . and use least action principle (ITUj) pQ 
d, (if ~ [U, [L», L v )]j + \mI{U - U- 1 ) = 0. 
In general case (nonstatic), energy of SU(2) Skyrme model is 



E = I d 3 x T 00 , 



where T 00 is energy-momentum tensor. Explicitly, 



E 



d x Tr 



F 2 1 
_ l~6~^ a ^ a ~ 32a 2 ^ c -^ a ' ^ c 



+^L L + -L^[L ,L a }[L ,L a } - ^M^U- 1 + U - 21) - ^L L C 

oCb 

L s tatic -L ro t a tj on 



where 



and 



E, 



static 



- / d A x Tr 



rp2 | p2 

— L 2 + \L a , L c } 2 + — MlOJ- 1 + U-2I) 

16 a 32a 2 1 ' J 16 nK ' 



Ei 



rotation 



- / dTx Tr 



16 



16a 2 



(10) 



'iV. 



(12) 



(13) 



(14) 
(15) 

(16) 



Solitonic properties of SU(2) Skyrme model for static energy in equation f fl5|) is studied 
by scaling the spatial coordinates 

x — > 2x/aF 
and express energy in F/Aa, i.e. by taking 

(F/Aa) = (1/127T 2 ). 

In this unit, equation (fTS]) becomes 

Ll + i ([La, L c f + ml(U-' + U - 21))] (17) 



where 



= 2M w /aF. 

Euler-Lagrange equation (fl2j) in static case becomes 



d a \ L a — — [L c , [L a , L c 



ni 



-{U -u- 1 ) = 0. 



(19) 



5 Scale Stability 



Let us look at scale transformation below 

x — > Xx. 

We find that L a currents, by scale transformation fl20l) . transform into 



L a {x) U-\\x)^^- = \L a {\x). 



(20) 



(21) 



The effect of scale transformation to static energy ffTTl) . by ignoring pion mass term (because 
pion mass is small), is 

E Mstatic = J E ° + XE Sky 

where E a is static chiral energy term and E^ y is Skyrme energy term. 
From equation (T2"2"|) . we obtain 



(22) 



dE [A] 



A=l 



A ; 



;E a + E. 



Sky 



+ E. 



Sky, 



A=l 



and 



d 2 £ [A] 



dX 1 



= E a \ x=1 = 2E a . 

A=l A 



The requirement for extremum condition is 

d£ [A] 
rfA 



0. 



(23) 



(24) 



(25) 
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We apply extremum condition, fl25|) . to equation ( 123]) . we obtain 

E a = E Sky , (26) 

which it shows 

E a > 0. (27) 

So that equation (EMI) fulfills condition 

d 2 £ [A] , . 

-^U > 0. (28) 

Equation (1251) is minimum stable condition which implies that static energy /I2ty) is stable 
against scale perturbation. 

6 Topological Charge 

Static energy of SU(2) Skyrme model can be expressed, by ignoring pion mass, as 

Static = ^ f d 3 x f-^j Tr (h a ± ^e abc [L b ,L c ]) =F J d 3 x e abc Ti [L a L b L c ] (29) 

where e abc is Levi-Civita symbol, e abc = #123 • At energy lower bound 

^static > B, (30) 

where 

B = --^2 y rf 3 xe a6c Tr (L a L b L c ). (31) 

5 integral, l[JJ\) . is independent of space-metric tensor, i.e. topological quantity which is 
known as topological charge of SU(2) Skyrme model. 

7 Static Skyrme Equation in Spherical Coordinate 

In spherical coordinate, (r, 6,<f)), static Skyrme equation has the following form [lj 

1 f 1 - „ 1 



- d r [ L r - - i — [L e , [L r , L e \\ + r2 s j n 2^[ L *' t Lr ' L ^ 



+ - \ |l[L r , [Le, L r ]] + ^T- e [L,, [LeM]}) (32) 



Its solution takes form 

U{r) = exp{iF a {r,e,<f>)a a ), (33) 

where 

F a (r,0,<j>)=g(r)n a (9,<l>) (34) 

is Skyrme ansatz, g(r) is profile function, n a , a = 1,2,3 is component of unit vector, n, in 
the internal space of SU(2), and a a is Pauli matrix. 
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8 Skyrmion Static Energy and Its Solution 

As a result of scaling spatial coordinates, by ignoring pion mass, static energy of SU(2) Skyrme 
model can be stated as [1] 



E. 



static 



12tt : 



r(47r) / r 2 dr 



-f- ) +^sin z g I 1+ I ^ 



dr 



dr 



+ — sin g 



From equation fl35|) . by using least action principle 

gE 'static = 0, 

we can derive Euler-Lagrange equation for profile function pQ 

d 2 g ^ 2 . 2 / dg^ 2 
dr 2 r 2 ® \ dr 



(35) 



(36) 



— sin 2g 



+ 1$ 



2\ 1 1 

- J - — sin 2g - — sin 2 g sin 2g = 0. (37) 
dr I \r I r z 



Equation (|37|) is second order of nonlinear differential equation. 

The solution of equation ( |37|) will be worked numerically with finite difference method. In 
order that static energy (|35i) has finite value at r = and r = oo, then profile function g{r) 
must fulfills boundary conditions 

g(0) = 7T, ff(oo) = 0. (38) 



9 Quantized Rotational Energy 

Skyrmion quantization is worked by involving time dependent of U (r) as 

U(r) -»- U(r,t) = A(t)U(r)A(t)\ 

where 

e SU(2)- mteiaah 
AA ] = A ] A = I 

A is time dependent unitary matrix. 

By using equation (1391) . rotational energy can be stated as [I], [12] 



E, 



rotation 



oo oo 

f J 2-2 27T /■, 2 . 2 

} dr r sin g H I dr r sin (7 

6 ./ 3a 2 / 



^\ 1-2 

— + — sin g 

dr I r l 



Tr 



(39) 



(40) 



R -m 



- I Tr ft 2 

2 



where 



with f2 is angular velocity matrix of Skyrmion, and 



1 = 2 



ttF 2 
~6 



2 poo 



dr r 2 sin 2 g + 



2tt 
3^2 



j 2-2 

dr r sin g 



is Skyrmion moment of inertia. 



dg\ 1.2 

— + -r sin g 

dr J r z 



(41) 
(42) 

(43) 
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10 Finite Difference Method 

Assume that g(x) is a continuous function of one variable. The value of this function is given 
only for discrete value, and equidistant x [20J: 



9k = g(x k ), (44) 

where x k = x + kAx. Let us define quantity 

Ag k = g k+1 - g k (45) 

which is called forward difference at point x k . Apply ( 145]) for higher order of forward difference, 
we get 

A 2 g k = Ag k+ i - Ag k = g k+2 - 2g k+1 + g k , (46) 
A 3 g k = A 2 g k+1 - A 2 g k = g k+3 - 3g k+2 + 3g k+1 - g k (47) 
and so on. General form of forward difference is: 

& r g k = z:= (-iy Q g k+r - t . (48) 

Backward difference is defined as 

V# fc = g k - g k -i (49) 
and higher order of backward difference is 

V 2 # fc = Vg k - V£f fc _i = g k - 1g k -\ + 9k-2 (50) 
so on. General form of backward difference is 

V0 fc = S[ =0 (-1) J Q ^_ r+J . (51) 

Following, it is defined central difference which is symmetric with x k . 

$9k = 9k+\/2 ~ 9k-i/2 (52) 
and higher order of central difference is 

S 2 9k = 9k+i ~ %gk + 9k-\ (53) 
so on. General form of central difference is 

5 r 9k = S[ =0 (-l) i Q g k - r/2+i . (54) 

If central difference is observed, k + 1/2 index is out of the value which is served by discrete 
value. In order to obtain the integer value, mean value operator, /i, is used to operate with 
g(x) function as central mean 

[dk+i/2 + gk-1/2] (55) 
9 



and higher order of central mean is 



H*9k = ^ [n9k+i/2 + Mk-1/2] 
so on. Central mean of central difference is defined as 



(56) 



(57) 



11 Numerical Solution of Skyrme Equation 

Skyrme equation is stated in finite difference form, i.e. by substituting equation (|49|) . (153 
into Skyrme equation ( 137]) then we obtain 



$ 2 9k 
Ar 2 



1 + — sin 2 g k 



Y9k 
Ar 



— sin 2g k 



+ 



Ar 



_ . 1 1 2 

- — -sm2g k — j sin g k sm 2g k = 0. 



By reposisition of equation (pgjh we get 

.._ ; ^ 2 



(^) 2 [^sin 2(/ fc ]-(^) (f) + ^sin 2^ fc + £ siri 2 g k sin 2<?* 



[1 + £ sin 2 



Ar 2 



-(Vg fc ) 2 sin 2g fc - (Vg fc )2(Ar)r + sin 2g k Ar 2 + ^ sin 2 g k sin 2^Ar 2 

r 2 + 2 sin 2 



(5? 



(59) 



(V^) 2 sin 2g k - (Vg k )2(Ar)r + sin 2g fc Ar 2 + ^ sin 2 g k sin 2g fc A? 

r 2 + 2 sin 2 g k 



The problems are initial values which are needed to be defined, i.e. 

g(r = 0), g(r = dr). 
Boundary conditions, which are given by Skyrme, are 

g{r = 0) = 7T, g(r — > 00) = 0. 



(60) 



(61) 



(62) 



Apply the boundary conditions (|62|) for initial values (IBTj) . It is suitable for g(r = 0), and 
g(r = dr) must be chosen, so that it is suitable for g(r — > 00) = 0. This method is known as 
shooting method [J3]. The ploting result, as shown in Figure 1: Profile function 4>(r), can be 
compared with trial and error method, i.e. 



g(x) = 4 arctan e r . 
In order to calculate static energy ( 135]) . the following integral is used 



(63) 
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4tan- 1 fe-^ 




Figure 1: Profile function <p(r) 



E{g) = ^ / r/r 



(r 2 g' 2 + 2 sin 2 g) + sin 2 g \2g' 2 + 



a , sm 9 



(64) 



where g' = dg/dr. Equation (IMj) . is solved numerically with trapezoid method and it gives: 

23.2154 Tr/^/e. 

Let us calculate numerically Static (|35|) . using numerical value of p(r) 



1 



static 



r(47r) 



23.2 



12vr 2V ~ "' V 2 

In the form of F and a, using definition (1/127T 2 ) = (F/Aa), we obtain 

vrF /23.2^ 



(65) 



E. 



static 



(66) 



a v 2 

The numerical value of g>(r) can be used to calculate 

(1) nucleon {m N ) and delta static masses (m^); 

(2) Skyrmion moment of inertia, I. 

We get benefit from numerical value of profile function, g(r), to calculate Skyrmion mo- 
ment of inertia. It gives 

1 /447\ 

' - 4 W) ■ < 67 > 

In quantum mechanics, angular momentum J is quantized as 

J 2 = j(j + l)h 2 (68) 
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where j = 0, §, 1,|,2,... 

2tt' 

/i is Planck constant. Here, we use using natural unit, h = 1 = c. 
Rotational energy is 

_ J 2 _ j(j + 1) 

-^rotation 

12 Hadron Mass 

Skyrmion is solution of Skyrme equation which has finite energy. Based on Einstein 
mass formula, Skyrmion mass is 

_ E static + E rotatlon _ 

m — o ~~ M static + ^rotation- 

C 2 

From equation (169]) and (!70|) . we obtain 

m = M stat i C H — — . 

Wess-Zumino quantization condition [3], requires: 

. _ 1 3 5 
J ~ 2' 2' 2' " ' 

It means that Skyrmion is fermion. 

In case of nucleon with spin |, j = |, equation ([71]) gives 

= M static + 2 ^ 2 ^|~ - 

By substituting the values of F — 123 MeV and a = 4.95 into (1ST)]) and (IBTj) . give 

M statJC = 905.08 MeV 

and 



I ^ 7.49 x 10~ 3 . 
Substitute ([73]) and ([74]) into ([72]) gives the result 

m N = 955.15 MeV. 

In case of delta particle with spin |, j = |, equation (17T]) gives 

+ 1) 

™a = M statlc + 2 2 2J ■ 

By substituting the values of M siaiic and / above, it gives 

m A ^ 1155.41 MeV. 
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13 Discussions 



The calculation of nucleon energy-mass gives the value about 955.15 MeV. The experimental 
value of nucleon energy-mass is about 939 MeV. The different value of nucleon energy-mass 
based on calculation and experiment is about 16.15 MeV. The calculation of delta energy-mass 
gives the value about 1155.41 MeV. The experimental value of delta energy-mass is about 
1232 MeV [3]. The different value of delta energy-mass based on calculation and experiment 
is about 76.59 MeV. 

The different values among calculation and experiment are becaused of SU(2) Skyrme 
model is very simple model, i.e. it consists of two flavours, (u, d) or (u,s) or (d,s). It is 
necessary keeping in mind related with these different values [T5] : 

(1) Meson field Lagrangian only includes pseudoscalar field. Other low mass meson (vector) 
should be included. 

(2) In nature, there are three flavour families than two "light" flavour families which must 
be included in more realistic formula. 

(3) The effects of chiral and flavour symmetry breaking aren't calculated yet. 

(4) Nc (number of colour) correction of nucleon mass isn't included yet. 
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15 Appendix: Flow Chart of Skyrme Profile Function 
Calculation and Its Integration 

clear; 
elf ; 

okmainloop=0 ; 
N=l; 

dr=.001; 
$%$dr=.01; 
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rf=50; 
°/orf=8; 
rf=50; 
rf=80; 

nnl=length( [0 dr:dr:rf rf +dr] ) ; 

g(nnl)=0; 

n0=l; 

n=2; 

V V V °/ V °/ V °/ °/ 

/o/o/o/o/o/o/o/o/o 

v=l+.0041 
vO=v; 

g(nO)=pi;°/„+0.01; 
g(n)=g(n0)-v*dr;°/ 0.02; 

if okmainloop 

% Looping utama program 

for r=[dr:dr:rf] 

n=n0+l; 

nl=n+l; 

gg=g(n) ; 
bg=g(n)-g(n0) ; 
bg2=bg*bg; 

sin2g=sin(gg)*sin(gg) ; 
sing2=sin(2*gg) ; 
r2=r*r ; 
dr2=dr*dr; 

% x=eFr (catto) 

g(nl)=(-2*r*bg*dr-sing2*(4*bg2-4*sin2g*dr2/r2-dr2))/(r2+8*sin2g)+2*g(n)-g(n0) ; 
n0=n0+l; 

°/ if (g(nl)<=0) break; end; 

end 

% akhir looping utama 

save simpan.mat 
else 

load simpan.mat 
end 

rr=[0 dr:dr:rf rf +dr] ; 
nn=nl ; 
gg=g(l:nl) ; 
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°/„subplot(2,l,l) ; 
subplot (1,1,1) ; 

plot(rr(l:nn-l),gg(l:nn-l),'b-'); 

rr=[0 dr:dr:rf rf+dr] ; 
°/ rre=pi*exp(-rr) ; 
rre=4*atan(exp(-rr) ) ; 

hold on; 

plot(rr(l :nn-l) ,rre(l :nn-l) , 'r- ') 
grid on 

/o /o /o /o /, /o /,°/o negatif change to zero 
°/„g=g.*(g>0); 

% Statical calculation 

IIEUc=0; 

IIaUc=0; 

IIbUc=0; 

IICUc=0; 

IIrUc=0: 



jumlah=0; 

dgdr2=v*v; % from initial value 

gg=g(D ; 

sin2g=sin(gg)*sin(gg) ; 

for n=[l:nnl-l] 
dg=g(n+l)-g(n) ; 
°/ dgl=g(n+2)-g(n) ; 
if n==l 

dgdr2=(dg*dg)/(dr*dr) ; 

else 

bg=g(n)-g(n-l) ; 
Dg=(g(n+l)-g(n-l))*.5; 
dgdr2=(dg*dg)/(dr*dr) ; 
°/ dgdr2=(dg*bg)/(dr*dr) ; 
°/ dgdr2=(Dg*Dg)/(dr*dr) ; 

end 

gg=g(n) ; 
r=n*dr ; 



sin2g=sin(gg)*sin(gg) ; 

j uml ah= j uml ah+dr * gg ; 
r2=r*r ; 

IIEUc=IIEUc+ dr*( r2*dgdr2 + sin2g* (2+8*dgdr2+4*sin2g/ (r*r) ) ); % catto 
IIaUc=IIaUc+ dr*r2*sin2g* ( 1 + 4*(dgdr2 + sin2g/(r*r)) ); °/„ catto 
IIbUc=IIbUc+ dr*r2*(l-cos(gg))*(l + dgdr2 + 2*sin2g/(r*r) ) ; 
IICUc=IICUc+ dr*r2* (1-cos (gg) ) ; 
IIrUc=IIrUc+ dr*r2*sin2g* (bg/dr) ; 
if (g(n+l)<0) break; end 

end 

IEUc=IIEUc*pi*4/8; 
IaUc=IIaUc*pi*2/3; 
IbUc=IIbUc*pi*l/2; 
ICUc=IICUc*pi*l/2; 
IrUc=IIrUc*2/pi; 

save II. mat II* dr vO 
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